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Abstract. Dynamical systems methods are used to study the evolution of the Bianchi I model with a scalar field. We show 
that inclusion of the non-minimal coupling term between the scalar field and the curvature changes evolution of the model 
compared with the minimally coupled case. In the model with the non-minimally coupled scalar field there is a new type 
of singularity dominated by the non-minimal coupling term. We examine the impact of the non-minimal coupling on the 
anisotropy evolution and demonstrate the existence of its minimal value in the generic case. 
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INTRODUCTION 

In the simplest cosmological model with a scalar field 
one can introduce a term modelling its coupling with the 
curvature. Such a term appears naturally if we understand 
the General Relativity as an effective field theory of 
gravity [1]. To answer the question how stable effects of 
the non-minimal coupling are, we consider the simplest 
anisotropic space of Bianchi I type which represents the 
simplest flat model with anisotropy. 

In the model under consideration we assume a uni- 
verse filled with a non-minimally coupled scalar field 
with the curvature and an unknown coupling constant ^ . 
The action assumes the following form 




- I |d V^g {eV> V a + e$R<j> 2 + 2tff»} 

(1) 

where K 2 — %TtG, £ = ±1 corresponds to a canonical and 
phantom scalar field, respectively, the metric signature 
is (—,+,+,+) and U(<j>) is the scalar field potential 
function assumed in a polynomial form. 

BIANCHI I METRIC 

We assume the following form of the metric 

ds 2 = -dt 2 + a 2 (t)(g 2 (t)dx 2 +g 2 (t)dy 2 +g 2 (t)dz 2 ) , 

(2) 

together with the condition 

glg2g3 = 1 • (3) 



The energy conservation condition we obtain from 
variation of the action (1) with respect to the metric 
components 

3H 2 -T.= K 2 (e-<j> 2 + U((l>) 

V 2 (4) 

+ ££(3// 2 -£)</> 2 + e3£H(0 2 )-J 
where the anisotropy is measured by 




the Hubble function is given by H = - and the space 
volume is V = a 3 . 

The acceleration equation is 

H = -2H 2 -^T.+ 

k 2 -e{l-6^ 2 + 4U{(p)~6^(j)U'((j)) (6) 
+ ~6 l-e£(l-6£)fc 2 </> 2 ' 

From space-space Einstein's equations we have addi- 
tional equations of motion for anisotropy 

(1 - e^K 2 <j>)q, = el; K 2 qi ^ 2 )~ 3H qi (l - k 2 <$> 2 ) . 

(7) 

When the non-minimal coupling is present one can 
write the right hand sides of Einstein field equations 
in several possible inequivalent ways [2]. In the case 
adopted here the energy momentum tensor of the scalar 
field is covariantly conserved, which may not be true 



for different approaches. For example, one can redefine which removes singularities of the system at z 2 = i 

.. . • . ..-9 . -9 _f i9 . . " ' tDI = 



the gravitational constant = K~ 2 — £<i;0 2 making it 
time dependent, then the effective gravitational constant 
can diverge for a critical value of the scalar field e = 
± 1 / \Jek 2 ^ and the model is unstable there with respect 
to arbitrary small anisotropy perturbations which become 
infinite there [3]. 

In what follows we introduce the energetic (expansion 
normalised) variables 



_ K (j) _ K VZTq 



zs ^' ^ = W 2 ' 



and 1/(0) — > f{z), then the energy conservation condi- 
tion can be expressed as 



1 =y 2 f(z) + e{l-6$)x 2 + e6${x + z) 2 + 
+ (l-e6^ 2 )i2 r , 

and the acceleration equation 

H 



(8) 



H 2 



e(l-6t;)x 2 -y 2 (2f(z)~3t;zf(z)) 



(9) 



l-£6^(l-6^)z 2 
The equations of motion for the anisotropy are 

(l-e6^ 2 )i^ = el2^ z -3(l-e6^ 2 ) (10) 
h qi 

Note that for the minimally coupled scalar field £, = 
equations (7) and (10) can be directly integrated and the 
evolution of the anisotropy does not depend whether the 
universe is filled with a canonical or phantom scalar field. 
In general case of the non-minimal coupling this is not 
possible because these equations depend on the phase 
space variables. 

The dynamical system is in the following form 



dx 
dlna 



--x(\-£l z )-e l -y 2 f{z)- 
-(x + 6§z)(^+2 + Q r ) 



H 



dlna y H 2 ' 
dz 



(11) 



dlna 



where using equations (8) and (9) we eliminate Q.^ and 
H/H 2 . One can see that right hand sides of the dynamical 
system (11) are rational functions of their arguments. 
In order to use standard dynamical systems methods we 
need to make the following time transformation 



(l-£6<?z 2 )(l-£6<?(l-6<?)z 2 ) d 



dlna dr 



(12) 



and z 



e6£(l-6£)- 



In this work we concentrate our attention at the fol- 
lowing critical point 



(x* = -6t;z*, f = 0, (z*) 2 



1 



e6S(l-6S) 



which for the isotropic case corresponds to the finite 
scale factor singularity [4]. The existence of this critical 
point depends on the value of the coupling constant E, 
only. For the canonical scalar field e = +1 it exists for 
< % < g while for the phantom scalar field e = — 1 for 
£<0or£<I. 

Eigenvalues of the linearization matrix at this critical 
point are 



36<g 2 
'1-6$' 



72<g 2 

i-6<r 



Stability, in time T, of this critical point depends on the 
value of the coupling constant | . For | < g the eigenval- 
ues are negative indicating a stable critical point, while 
for £ > 1 the eigenvalues are positive which indicates an 
unstable critical point. 

The solutions of the linearised system in the vicinity 
of the critical point under investigation are 



x(t) = x* + (Ax + 2(1 -3^)Az)exp(AiT)- 

-2(l-3^)Azexp(A 3 T) 
>'(t) =y*+Ayexp(A2T), 
z(t) =z*+Azexp(A 3 r). 



(13) 



where Ax = - x* , Ay = - y* and Az = - z* 
are displacements of the initial conditions with respect 
to coordinates of the critical point. 

From (8) we see that at this critical point = 0, and 
using the linearised solutions we have 



Q e (t) - £2(1 -6^) 2 z*Azexp(A 3 T) 



(14) 



On the other hand the anisotropy parameter E = 3H 2 £lz 
and in the energy phase space variables ^Tj^£ = ^r- 
Using linearised solutions in the vicinity of this critical 
point we have 



1 



K 2 U Q 



Qe(t) 
(y(T)) 2 



2 * 



e2(l-6^) 2 z 



Az 
(Ay) 2 



const. 



(15) 

As long as the linear solutions are considered, one can 
see, that the anisotropy function is strictly constant and 
the value of anisotropy in the vicinity of the critical point 
representing the finite scale factor singularity depends 
on the value of the non-minimal coupling constant E, 
and initial conditions in the linear approximation. Of 
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FIGURE 1. 

calculated from (8) and using 3H Z 



Evolution of the anisotropy function £ = -Mj- 

j2 _ eua. The initial con . 

y 

ditions taken in the vicinity of the critical point under inves- 
tigations, and U((j>) °c (j) 2 , e = +1, £, = 5. We observe that 
the anisotropy is constant at the critical point then approaches 
a minimal value and next grows as the universe expands. A 
dashed vertical line denotes value of the scale factor of the 
singularity. 



course one can feel bad about the square value of Ay 
in the denominator which should vanish in the linear 
approximation, but in the dynamical system (11) the 
phase space variable y appears only in a polynomial form 
of the second degree and one can make the following 
change of variables y—y 2 and then the same results hold. 

In Figure 1 we plotted evolution of the anisotropy. At 
the first plot we present evolution of £ with respect to 
time T and the critical point representing a finite scale 
factor singularity is approached at T — > °°. At the second 
plot we present evolution of the anisotropy with respect 



to the scale factor a/ a ^ rescaled to the value of the scale 
factor taken at the initial conditions. 



KASNER METRIC 



In this section in (2) we assume 



gl=t 
and the metric is 

,2 



82 = t' 



S3 



- tP3 



ds 2 = -df 2 - 



■a 



(16) 



(f ) (t 2pi dx 2 + t 2pi -dy 2 + t 2p3 dz 
where from (3) we have constraint equation 
Pi +P2 + P3 =0. 



The assumed form of the metric represents a special case 
of the general Bianchi I metric which can be recast into 
the standard Kasner solution for an empty space. 
Now the anisotropy function is 



■pI 



2^1 ,° L 



(17) 



and the energy conservation condition expressed in the 
energetic variables is 



1 =y 2 f(z) + e(l-6i;)x 2 + e6i;(x + z) 2 - 
+ {\-e^z 2 )o 2 h 2 7 

where we denned 

h : 



(18) 



1 

7Tt 



From (10) we obtain 

(1 - e6%z 2 )h = 3(1 - e6$z 2 ) - ell^xz. 
The dynamical system is in the following form 



(19) 



dx 
dlna 



dy 

dlna 

dz 
dlna 



1 



-x(l -a 2 h l ) -e-y 2 f (z)- 
f H 

-(x + 6£z)[jZ2 +2 + a 2 h 2 
H 



(20) 



where using equations (18) and (9) (after substitution 
ill = 2 h 2 ) we eliminate terms containing a 2 h 2 and 
H/H 2 . The dynamical system describing evolution in 
the Kasner metric is exactly in the same form as in the 
case of the general Bianchi I model. One exception is 
the existence of the additional constraint equation (19). 
Using this equation we find that at the critical point under 



where /; is given by (19) and using the linearised solu- 
tions obtained in the previous section we find 



0.008 
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a 2 = e2(l -6|)VAzexp(A 3 T). (22) 

In Figure 2 we plotted evolution of the anisotropy pa- 
rameter a 2 for the Kasner metric. In evolution in time 
T (first plot) critical point under investigations is ap- 
proached as T — > °°. At the second plot we presented 
evolution of O 2 with respect to the scale factor a /any 
Both plots show that during the contraction of universe 
the anisotropy for the Kasner metric decreases and van- 
ishes at the finite scale factor singularity. 

CONCLUSIONS 

Dynamics of the simple anisotropic cosmological model 
with non-minimally coupled scalar field was investigated 
near the singularity. The critical point representing the fi- 
nite scale factor singularity exists only for non-minimal 
% 7^ and non-conformal | ^ \ values of the coupling 
constant. We have shown that for the most general form 
of the metric anisotropy of the universe reaches mini- 
mal value and then model starts to anisotropies. More- 
over in the finite scale factor singularity the value of 
the anisotropy tends to a constant value which depends 
on the non-minimal coupling. On the other hand, if we 
consider only small anisotropies, the universe isotropises 
during the contracting phase and reaches an isotropic flat 
cosmological model in the finite scale factor singularity. 
As the universe expands small anisotropies grow indicat- 
ing instability with respect to arbitrary small anisotropy 
perturbations. 



FIGURE 2. Evolution of the anisotropy parameter a 2 calcu- 
lated from (21), and U(<j>) <* (j) 2 , £ = +1, % = \. In the Kas- 
ner metric anisotropy vanishes at the critical point represent- 
ing the finite scale factor singularity. As the universe expands 
anisotropy grows indicating that this critical point is unstable 
with respect to anisotropy perturbations. A dashed vertical line 
denotes value of the scale factor of the singularity. 



investigations h* = 1 this indicates that while a sample 
trajectory approaches a critical point and H — »• oo then 

The energy conservation condition (18) can be rewrit- 
ten into the following form 

l-y 2 f(z)-e(l-6^x 2 -e6^x + z) 2 .... 

(l-e6<^ 2 )/z 2 ( ' 
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